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We derive equations for the mean entropy and the mean internal energy in the low- 
Mach-number temperature stratified turbulence (i.e., for turbulent convection or stably 
stratified turbulence), and show that turbulent flux of entropy is given by F s = pus , 
where p is the mean fluid density, s are fluctuations of entropy and overbars denote aver¬ 
aging over an ensemble of turbulent velocity field, u. We demonstrate that the turbulent 
flux of entropy is different from the turbulent convective flux, F c = Tpus, of the fluid 
internal energy, where T is the mean fluid temperature. This turbulent convective flux 
is well-known in the astrophysical and geophysical literature, and it cannot be used as a 
turbulent flux in the equation for the mean entropy. This result is exact for low-Mach- 
number temperature stratified turbulence and is independent of the model used. We also 
derive equations for the velocity-entropy correlation, us, in the limits of small and large 
Peclet numbers, using the quasi-linear approach and the spectral r approximation, re¬ 
spectively. This study is important in view of different applications to the astrophysical 
and geophysical temperature stratified turbulence. 


1. Introduction 

Temperature stratified turbulence (e.g. , turbulent convection or stably stratified turbu¬ 
lence) plays a crucial role in astrophysics (Shakura et oil 1 978t Peebles 1980; Zeldovich et al 


1983, 1990; Ruziriaikin ef ai.lll988l: IC'larke fc Carswelll2007 ) and geophy sics ( Monin fc Yagloml 


197.5c 


Zilitinkevichl ll99lT Canuto 20091: Zilitinkevich et al. 2008L 2013l l. The large-scale 


properties of the temperature stratified turbulence are determined in the framework of 
the mean-field approach in which all quantities are decomposed into the mean and fluc¬ 
tuating parts, where the fluctuating parts have zero mean values and overbars denote 
averaging over an ensemble of turbulent velocity fields. 

In the astrophysical and geophysical literature on the low-Mach-number temperature 
stratified turbulence two different formulae for the turbulent flux of entropy are used. The 
first formula coincides with the turbulent convec tive flux of internal energy, F r = Tp us 


so that the equation for the mean entropy is (IKi tchatinov & Mazur 2000; Bru n et al 


2004: Miesch et al. 2008: Jones et al. 200 


o pv is IIKit cliatmov &; IVlaz uill/UUUt Inrun et a,L 

Jones fc Kuzanvai]20091: Kapvla etHIl 20121): 


^ + (U • V)s) + • (Tp us) = = [V • (KVT) + J] , (1.1) 


where p, T, S and U are the mean fluid density, temperature, specific entropy and 
mean velocity, respectively; u and s are the fluctuations of fluid velocity and entropy, 


f Email address for correspondence: gary@bgu.ac.il 

























































2 


I. Rogachevskii and N. Kleeorin 


respectively; J is the mean source and/or sink of the entropy (that also includes the 
viscous heating) and K is the coefficient of the molecular heat conductivity. The last 
term on the left-hand side of Eq. (O) corresponds to the turbulent flux of entropy. 

The other for m of the turbulent flux of entropy is F„ = pTZs, and the equ ation for the 
mean entropy is ( Braginsky fc Robertdll995 : Glatzmaier fc Roberts! 1996«M ): 


-p + (u • V)s) + V • {pus) = = [V • (. KVT) + J] . (1.2) 

Equations ED and (11.21) are essentially different. In particular, the last term on the 
left-hand sides of Eq. ED and ED are different. 

The goal of the present paper is to derive equations for the mean entropy and the 
mean internal energy which yield formulae for the turbulent flux of entropy and the 
turbulent flux of internal energy, and to clarify which equation for the mean entropy 
[Eq. (11.11) or ED] used in the temperature stratified turbulence, is correct. When the 
fluid temperature profile is not uniform, the above question is crucial. 


2. Turbulent convective flux of mean internal energy and turbulent 
flux of mean entropy 

In this Section we will derive equations for the mean entropy and the mean internal 
energy. We consider low-Mach-number temperature stratified fluid flows. 


2.1. Governing equations 


The budget equation for the instantaneous internal energy density E 
1959th 


c v T is (Landau & Lifshitz 


d{ P E) 

dt 


+ V-{pUW 


UP - KVT) 


Q, 


( 2 . 1 ) 


where U is the instantaneous velocity determined by the Navier-Stokes equation for fluid 
motion, p , T and P are the instantaneous density, temperature and pressure, respectively, 
which satisfy the equation of state for a perfect gas, K is the coefficient of the molecular 
heat conductivity, W = c p T = c v T + P/p = E + P/p is the instantaneous enthalpy, 
where c v and c p are the specific heats at the constant volume and pressure, and Q = 
-PV • U+aij(U)S7jUi, where a{U) = 2upS{U ), 5(17) = 5y = \{U^+Uj^-\SijV-U, 
v is the kinematic viscosity, and Sij is the Kronecker tenzor. The instantaneous density 
p is determined by the continuity equation: 


| + V.(,V)=0. 


( 2 . 2 ) 


The budget equation for the instantaneous kinetic energy density \pU 2 is 
195911 : 


( Landau fc Lifshitj 


dt 


(\pU 2 ) + V • [U (\ P U 2 + P) - Ua(U)] = -Q. 


(2.3) 


The sum of Eqs. ED and (12.31) yields the conserva tion law for the instanta neous total 
(kinetic plus internal) energy densities \pU 2 + pE ( Landau fc Lifshi~t3 1959fl : 


— (±pU 2 + pE) + V ■ [U {\pU 2 + pW) - Ua{U) - KVT] = 0. 


(2.4) 
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2.2. Turbulent convective flux and equation for mean internal energy 

Averaging Eq. HU) over the ensemble we obtain the budget equation for the mean 
internal energy density E = c v T: 


d(pE) 

dt 


+ V • (pU E + p uw — up — K'VT'j = Q , 


(2.5) 


where u are the velocity fluctuations, w = T s + p/p are the enthalpy fluctuations, s 
and p are the entropy and pressure fluctuations, respectively, and Q = — PV • U — 
pV • u + <7,j(t/)V jUi + &ij(u)V jUi. In the derivation of Eq. (12.51) we used the identity: 
W = E + P/p, and since we con sider a lo w-Mach -number turbulent flow, we took into 
account that |wp|/|p| <C |Us|/|S'| (jChassaing et al. 120021 ). The turbulent flux of enthalpy 
is 


_ —_ up 

uw = T us H-. 


( 2 . 6 ) 


Since W = E+P/p we obtain that pU E = U ( pW — P ). Substituting the latter equation 
and Eq. (12.61) into Eq. (12.51) . we obtain: 


+ V- (pUE + Tpus-KVT) =Q. (2.7) 

The mean fluid density p is determined by the continuity equation: 

% + V ■ (p U) = 0. (2.8) 

Equations m and (12.81) yield the following equation for the evolution of the mean 
internal energy: 

P + {U ■ V)E^j +V-(Tpus- KVT) = Q, (2.9) 

where the turbulent convective flux of the mean internal energy is: 

F c = T~pus. (2-10) 


2.3. Equation for the sum of mean and turbulent kinetic energies and conservation law 

for total mean energy 

Averaging Eq. (12.31) for the instantaneous kinetic energy density \pU 2 we obtain an 
equation for the sum of the mean and turbulent kinetic energies \~pU + ^ ~pu 2 : 


d_ 

dt 


(\pU~ + \pu 2 ^j + V • U [\pU 2 + P^j — Ucj{U) + u ( \pu 2 +p) — ua(u) 


= -Q. 

( 2 . 11 ) 

The sum of Eqs. m and (12.111) yields the conservation law for the total mean energy 

E to t = \fU 2 + Ppu 2 + pE: 

d 


— [pU 2 + ^pu 2 +pE^ + V- U (pU 2 + P + -pE') ~Ua(U) 


—uir(u ) + u (\pu 2 + p) + Tp us — K"VT 


= 0 . 


( 2 . 12 ) 


This equation contains the turbulent convective flux F c = Tpus. The conservation 
law (12.121) for the total mean energy Etot can be rewritten in terms of the mean, U W, 
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and turbulent, mu, fluxes of enthalpy [see Eq. my 

| ( W 2 + + pE)+v-\u (\-pu 


pW 


— ua(u) + u (^pu 2 ) + puw — K'VT 


= 0 . 


Ua(U) 


(2.13) 


2.4. Equation for mean entropy 


T he e volution equation for the instantaneous entropy S = c v ln(Pp 7 ) is (Landau & Lifshitz 
1959): 


d_ 

dt 


U V 


s = -f [V • (A'VT) + J] 


(2.14) 


where 7 = c p /c v is the ratio of specific heats and J is a source and/or sink of the entropy 
(that also includes the viscous heating). Multiplying the equation for the entropy (12.141) 
by the fluid density p, and the continuity equation m by the fluid entropy S, and add 
them, we obtain the following equation: 

+ V • (P US) = i [V • (A'VT) + J]. (2.15) 

The second term on the left-hand sides of Eqs. (12.141) and (12.151) . which contributes to 
the turbulent diffusion of the mean entropy, does not contain the temperature field. This 
is a reason why the turbulent flux of the mean entropy for stratified turbulence cannot 
contain the mean temperature. 

Averaging Eq. (12.151) over the ensemble we obtain the equation for the mean entropy, 
S: 


d(pS) 

dt 


+ V -(pUS- 


pus)= = [V-(KVT) + J} 


(2.16) 


In the derivation of Eq. (12.161) we have taken into account that for a low- Mach-n u mber 
turbulent flow: |ps| <C |p[ IS), |ps|/|p| -C |tts|/u r ms and |up[/|p| -C |Tts|/|S| (jChassaing et al. 
2002f) . To get the simplest form of the molecular diffusion term and the source term on 
the right-hand side of Eq. (12.161) we assumed that: (a) the temperature fluctuations, 0, 
are much smaller than the mean fluid temperature, T, i.e., |0| (b) in the frame¬ 

work of the mean-field theory there is a separation of scales, £ 0 -C Tt, where Lt is the 
characteristic scale of the mean temperature variation and £0 is the integral scale of tur¬ 
bulence (the random velocity field); (c) the coefficient of the molecular heat conductivity 
K is independent of the temperature fluctuations and (d) fluctuations of the source or 
sink of the entropy, J 1 , are independent of the temperature fluctuations. Equations (12.81) 
and (12.161) yield the following equation for the evolution of the mean entropy: 


P + (u • V)s) + V • (pus) = i [V • (KVT) + J] , (2.17) 

which coincides with Eq. m, and the turbulent flux of the mean entropy for stratified 
turbulence with non-uniform profiles of the mean fluid temperature and density is: 


F s = pus. (2-18) 

Other forms of the turbulent flux of entropy [see Eq. lO] used in the astrophysical and 
geophysical literature are incorrect. This is an exact statement for a low-Mach-number 
temperature stratified turbulence and is independent of the model. 
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3. The velocity-entropy correlation 

Now let us determine the velocity-entropy correlation, us. For simplicity we consider 
turbulent flows with a zero mean velocity, U = 0. Subtracting Eq. (12.1611 from Eq. (12.151) 
we obtain the equation for the entropy fluctuations: 

r) q 

—+A r- x v 2 s = l, (3.1) 

where M = (p) _ 1 V- [p(us — Its)] is the nonlinear term, y = K/cfp is the molecular 
diffusion coefficient of entropy and I = — (p) _1 V • (~pS u) = — (it ■ 'V)S is the source 
term. In the derivation of Eq. <33* we took into account the anelastic approximation, 
[V • (jsu) = 0], we assumed that fluctuations of the source or sink of the entropy, J', 
are very small, and we also assumed that the molecular diffusion term can be simplified 
as (pT) _1 V • [(A'/c p )V(Ts)] ~ yV 2 s. In the latter estimate we assumed that: (i) the 
temperature fluctuations are much smaller than the mean fluid temperature; (ii) £q Ft! 

(iii) the coefficient of the molecular heat conductivity is independent of the coordinates; 

(iv) for a low Mach numbers the entropy fluctuations are given by 


s = c p 



(i -t)p \ 

C sP ) 



(3.2) 


_ ^ j 2 __ 

where c s = (yP/p) ~ is the sound speed. Equation (13.21) follows from the equation of 
state for a perfect gas: P = (c p — c v )pT and the definition of entropy. Let us derive the 
equation for the velocity-entropy correlation us in two limiting cases for small and large 
Peclet number, where Pe = uq£o/x is the Peclet number and u o is the characteristic 
turbulent velocity in the integral scale of turbulence, £q. 


3.1. Small Peclet numbers 

In order t o study entropy fluctuations for small Peclet numbers we use a quasi-linear 
approach ( Moffatt 1978 : Krause fc Raedleill 198(1) . that for a given velocity field is valid 
only for small Peclet numbers (Pe <C 1). In the framework of this approximation we 
neglect the nonlinear term and keep the molecular diffusion term in Eq. <33*. We rewrite 
Eq. m in a Fourier space and solve this equation. The solution is: 


s(u>, k) = G x (w, fc)/(w, k), (3.3) 

where G x (uj, k) = (yfc 2 + iw) _1 , w is the frequency, and k is the wave vector. We apply 
a standard two-scale approach, whereby the non-instantaneous two-point second-order 
correlation function is written as follows: 


Ui(ti,x) s(t 2 ,y) 


Ui(wi,ki)s(w2, fe 2 ) exp^fei-a: + k 2 -y) 


+i(uiti + u} 2 t 2 )\ duii dui 2 dk\ dk 2 


J Fi(ui, k) exp[ifc-r + iuif] duidk, (3.4) 


where we use large scale variables: R = (x + y)/ 2, K = k\ + k 2 , t = (ti + t 2 )/2, 
12 = uq + ui 2 , and small scale variables: r = x — y, k = (fci — k 2 )/2 , t = ti — t 2 , 
ui = (uji — uj 2 )/2 , and 


Fi(ui, k) 


Ui(uii, ki) s(ui 2 , k 2 ) exp[iflt + iK-R] dQ dK. 


(3.5) 


Here cci = uj + £2/ 2, ui 2 = —u> + ST/2, k\ = k + K/ 2, and k 2 = —k + K /2 (see, e.g., 
Roberts fc Sowar J (11975 V). We assume here that there is a separation of scales, i.e., the 
maximum scale of random motions £q is much smaller than the characteristic scales of 
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inhomogeneities of the mean entropy and fluid density. Equations (13.3M3.5I) yield the 
velocity-entropy correlation us: 


UiS = J Ui(uj,k) I(—oj,—k)G x dw dk = —(V jS) j Ui(ui,k) Uj(—u},—k) G* dw dk. 


(3.6) 


We use the simple model for the second moments, ui(uj,k)uj(—ui,—k), of a random 
velocity field in a Fourier space for inhomogeneous, isotropic and incompressible flow: 


i(ui , k) —k) = 


E(k) $(w) 


8-7T k 2 


Sij 


ki kj i 

~F r + 2k ^ 


(fciVj kjVi 


u* 


(3.7) 


This model is obtained using the symmetry arguments ( Batcheloinil97ll h Here Sij is the 


Kronecker tensor, the energy spectrum function is E{k) = Ce k^ 1 (q — 1) ( k/ko)~ q for 
the range of wavenumbers ko < k < kd, the wavenumber ko = l/£o> the length £q is the 
maximum scale of random motions, the exponent 1 < q < 3, and the coefficient Ce = 
[1 — {ko/kd ) q ~ l ] -1 . We use the frequency function <f>(w) in the form of the Lorentz profile: 
<£>(w) = [7rr c (w 2 + t” 2 )] -1 , where t c is the correlation time of a random velocity field. 
This model for the frequency function corresponds to the following non-instantaneous 
correlation function: Ui(t)uj(t + r) oc exp(—r/r c ). 

We use Eqs. m and (ESP , and after integration in w-space and in fc-space in Eq. (1X61) 
we obtain the formula for the velocity-entropy correlation Ui s: 


UiS= -XT^iS , 

XT = C x uq £ 0 Pe, 


q~ 1 


1 - (VM 9+1 


° x 3(g + l) ’ 


(3.8) 

(3.9) 


where Xt is the coefficient of the turbulent diffusion of the mean entropy and «o = V 
is the characteristic velocity in the maximum scale of random motions. Here we used 
that I 0 = J $(w)G x (w, k) du = r c /( 1 + r c ??fc 2 ), and for small Peclet numbers Iq « 
(r/k 2 )^ 1 . The coefficient C x = 1/3 for a narrow range of the random velocity field in 
the wavenumbers, kd — ko <C kd, and C x = (q — l)/3 (q + 1) for a wide range in the 
wavenumbers, kd 3> &o- Contributions (which are proportional to V in Eq. (13. 71) ') to 
the velocity-entropy correlation uCs, after the integration over the angles in fc-space, 
vanish. However, the coefficient of the turbulent diffusion xt depends on the coordinates, 
due to the inhomogeneous turbulence. Equations (13.81) a nd (|3.9I) are in agreement with 


those obtained by means of dimensional arguments ([Batcheloretmi. 


Lagrangian-history direct-interaction approximation (jKraichnanl 1968). 


19591) and by the 


3.2. Large Peclet numbers 

In this subsection we derive a formula for the velocity-entropy correlation us using the 
spectral r approach that is valid for large Peclet and Reynolds numbers (Pe 1). Using 
Eq. (13.ID written in a Fourier space we derive equation for the instantaneous two-point 
second-order correlation functions Fi(t, k) = m(t , k) s(t, —k): 

^ = u 1 {t,k)I(t,-k) + MF^ III \k), (3.10) 

where MF^ III \k) = —[uiJV + ( dui/dt ) s — x^i V 2 s]fc are the third-order moment terms 
appearing due to the nonlinear terms which also include the molecular diffusion term. 

The equation for the second moment includes the first-order spatial differential opera¬ 
tors applied to the third-order moments. A problem arises regarding how to close the sys- 





























Turbulent fluxes of entropy and internal energy 7 

tern, i.e. , how t o express the third-orde r terms M.F ( JJJ ) through the lower moments F' 11 ' 
( Orszad 197f)l : Monin fc Yagloml 1975 : McComt J1 1 99()l i . We use the spectral r approxi¬ 
mation which postulates that the deviations of the third-moment terms, AiF^ 111 ^ (k), 
from the contributions to these terms by the background turbulence, A4F l ' III ' 0 \k), can 
be expressed through similar devi ations of the second moments, F( TI >(k) — F^ II,0 \k) 
( Orszad 19701 Pouauet et al. 1976h : 


1 


MF ( ~ III \k ) - MF^ IIIfi \k) = -— F^ u \k) - F^ n ’°\k) 

T r [k) L 


(3.11) 


where T r {k) is the scale-dependent relaxation time, which can be identified with the cor¬ 
relation time r(fc) of the turbulent velocity field for large Reynolds and Peclet numbers. 
The functions with the superscript (0) correspond to the background turbulence with 
a zero gradient of the mean entropy. Validation of the r approximation for different 
situations has been perfor med i n num erou s nu merical s imulations and analytical stud- 
ies (see e.g. the review by iBrandenburg &; Subramanianl (2005 1: and also discussions by 
Rogachcvskii V Kleeorin ( 2007 ): iRogachevskii et al. ( 2011 )). 

Note that the contributions of the terms with the superscript (0) vanish because when 
the gradient of the mean entropy is zero, the turbulent heat flux and the entropy fluc¬ 
tuations vanish. Consequently, Eq. (13.111) for MFj; III \k ) is reduced to MF^ III \k) = 
— Fi(k)/r(k ). We also assume that the characteristic time of variation of the second 
moment Fj(fc) is substantially larger than the correlation time r(fc) for all turbulence 
scales. Therefore, in a steady-state Eq. (13.101) yields the following formula for the velocity- 
entropy correlation: 

Fi = f r(k) (ui(t,k) I(t, —fc)) dk = —(V jS) f r(k) (ui(k) Uj(—k)) dk. (3.12) 


We use the following simple model for the second moments, Ui{k ) Uj(—k), of a turbulent 
velocity field in Fourier space for inhomogeneous, isotropic and incompressible flow for 
large Reynolds numbers: 


Ui(k) Uj(—k) = 


m 


Sin - 


hi kj 


[kiS7 j kj\7i ) 


(3.13) 


8t rfc 2 l 13 k 2 ' 2 k 2 

This model is obtained using symmetry arguments dBatchelor 1ll971 ). After integration 
in fc-space of Eq. (13.121) we arrive at an equation for the velocity-entropy correlation, 
(ui s): 


(ui s) = -XT^iS , xt = wo 4 / 3 , 


(3.14) 


where uq = V u 2 is the characteristic turbulent velocity. In the derivation of Eq. (13.141) 
we used the following expression for the turbulent correlation time: r(k) = 2 tq (k/ko) 1 ~ q , 
where to = 4/wo is the characteristic turbulent time. Contributions (which are propor¬ 
tional to V in Eq. (13.13D ) to the the velocity-entropy correlation ufs, after the integration 
over the angles in fc-space, vanish. However, the coefficient of turbulent diffusion xt de¬ 
pends on the coordinates, due to the inhomogeneous turbulence. Therefore, the formulae 
for the velocity-entropy correlation, (Ui s ), are similar for small and large Peclet num¬ 
bers, while the coefficients of turbulent diffusion of the mean entropy are different in 
these two limiting cases. Eq uatio n (13.141) is in agreement with that derived by means of 
the path integral approach ( Elnerin et all ll995l) . by dimensional arg ume nts a nd by the 
renormalization procedure used for large Peclet numbers ( Elnerin et al. 199fil ). 
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4. Conclusions 

In the present study we have demonstrated that for a low-Mach-number compressible 
fluid flow the turbulent flux of entropy, F s = ~pus = —J>xt^S, is different from the 
turbulent convective flux of the mean internal energy, F c = T~pus = —TJ>xt^S. As 
follows from the analysis performed in Section 3, the coefficient of turbulent diffusion 
of entropy Xt depends on the Peclet number. For small Peclet numbers, applying the 
quasi-linear approach for an isotropic and inhomogeneous background random velocity 
field we obtain the following coefficient of turbulent diffusion of the mean entropy: Xt = 
C^Petto^O: where the constant C x depends on the energy spectrum of the random velocity 
field. For large Peclet and Reynolds numbers, applying the spectral r approximation we 
get the following coefficient of turbulent diffusion of the mean entropy: Xt = uo^o/3- 


We are indebted to A. Brandenburg, E. Dorrny, P. J. Kapyla, M. Rheinhardt and 
P. H. Roberts for stimulating discussions. This work was supported in part by the Re¬ 
search Council of Norway under the FRINATEK (grant No. 231444) and the Academy 
of Finland (grant No. 280700). 
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